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Abstract 

Top quark production in the process e“*“e“ -A- tt at a, future linear electron positron collider with 
polarised beams is a powerful tool to determine indirectly the scale of new physics. The presented study, 
based on a detailed simulation of the ILD detector concept, assumes a centre-of-mass energy of ^/s = 
500 GeV and a luminosity of £ = 500 fb“* equally shared between the incoming beam polarisations of 
Ve-, 'Pe+ = ±0.8, =f 0.3. Events are selected in which the top pair decays semi-leptonically and the cross 
sections and the forward-backward asymmetries are determined. Based on these results, the vector, 
axial vector and tensorial CP conserving couplings are extracted separately for the photon and the Z'^ 
component. With the expected precision, a large number of models in which the top quark acts as a 
messenger to new physics can be distinguished with many standard deviations. This will dramatically 
improve expectations from e.g. the LHC for electro-weak couplings of the top quark. 


1. Introduction 

The main goal of current and future machines at the energy frontier is to understand the nature of 
electro-weak symmetry breaking. This symmetry breaking can be generated by the existence of a new 
strong sector, inspired by QCD, that may manifest itself at energies of around 1 TeV. In all realisations 
of the new strong sector, as for example Randall-Sundrum models [1] or compositeness models [2], the 
strength of the coupling to this new sector of the Standard Model fields are supposed to increase with 
their mass. For this and other reasons, the heavy top quark or t quark with a mass of approximately 

= 173 GeV [3] is expected to be a window to any new physics at the TeV energy scale. New physics 
can modify the electro-weak ttX vertex described in the Standard Model by Vector and Axial vector 
couplings V and A to the vector bosons A = 7 , At the International Linear Collider, ILC [4], 
that will collide electron and positrons at a centre-of-mass energy of 500 GeV, t quark electro-weak 
couplings can be measured at the % level. 

In contrast to the situation at hadron colliders, the leading-order pair production process e“''e“ —>■ tt 
goes directly through the ttZ'^ and ttq vertices. There is no concurrent QCD production of t quark 
pairs, which increases greatly the potential for a clean measurement. A parametrisation of the ttX 
vertex valid to all orders of perturbation theory may be written as * : 
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with e being the electrical charge of the electron, = {q-\- q)^ being the squared four-momentum of 
the exchanged boson and q and q being the four-vectors of the t and t quark, respectively. Further, 7 ^ 
are the Dirac matrices leading to vector currents of fermions and 75 is the Dirac matrix allowing to 
introduce an axial-vector current into the theory. Finally, | — Ivlti) allows for describing 

the scattering of a particle with spin 1/2 and a given magnetic moment. 

Within the Standard Model the Fi have the following values at tree level: 


_ 2 p-i,SM _ „ t^Z,SM 



riZ.SM 

^lA 




( 2 ) 


while all the F 2 are zero. In Eq. 2 Sw and Cw are the sine and the cosine of the Weinberg angle Ow- The 
scale dependence of the form factors is a consequence of higher order corrections. The corrections of 
the vector currents lead to the anomalous electro-magnetic and weak-magnetic moments represented 
by F^ that correct the gyromagnetic ratio gt of the t quark. Typical values for these corrections are 
in the range 0(10“^ — 10“^) [G]. Corrections to the axial-vector current result in the Form Factors 
that are related to the dipole moment df = (e/ 2 mt)F^( 0 ) that in turn violates the combined 
Charge and Parity symmetry CP. Otherwise said, all couplings but F^{k'^) conserve CP. 

The Form Factors Ffy ^ are related to couplings of t quarks with left and right-handed helicity to 
the Z°: 


gf = F^y - 5 I = F^y + (3) 

Trivially, the same equations apply correspondingly to the photon couplings g/ 

In this paper the precision of CP conserving form factors and couplings as introduced above 
will be derived by means of a full simulation study of the reaction e+e” -A ti & centre-of-mass 
energy of ^/s = 500 GeV with 80% polarised electron beams and 30% polarised positron beams using 
experimentally well defined observables. Special emphasis will be put on the selection efficiency and 
the polar angle of the final state t quarks. Both experimental quantities are suited to monitor carefully 
experimental systematics that may occur in the extraction of form factors and couplings. The results 
presented in the following are based on the studies described in detail in Refs. [7, 8]. 


2. Top quark production at the ILC 

The tree level diagram for pair production of t quarks at the ILC is presented in Figure la. The 
decay of the top quarks proceeds predominantly through t -A W^b. The subsequent decays of the 
bosons to a charged lepton and a neutrino or a quark-anti-quark pair lead to a six-fermion final 
state. The study presented in this article focuses on the ’lepton+jets’ final state l^vhhq'q representing 
a branching fraction of about 43.4% on all ti pair decays. 

Several other Standard Model processes give rise to the same final state. The most important 
source is single t quark production through the process e+e“ — WW* -A Wtb -A l^ubbq'q. One of 
the diagrams contributing to this process is presented in Figure lb. Another relevant source is the 
Z^W'^W~ production. Due to the coupling of initial state electrons or positrons to W bosons both 
sources contribute nearly exclusively in a configuration with left-handed polarised electron beams and 
right-handed polarised positron beams. 

In that case single t quark and Z'^W^W~ boson production can yield a total production rate of 
up to 10% of that of the pair production diagram of Fig. la. Experimentally, Z^W~^W~ production 
can be distinguished rather efficiently from ti pair production, but a clean separation of final states 
with a single t quark seems impossible. A realistic experimental strategy must therefore consider the 
W~^bW~b inclusively [9]. 
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(a) tt pair production 

Figure 1: Diagrams that contribute to the e~^e~ —>■ 
(b) single t quark production. 



(b) Single t quark production 
q production: (a) Tree level tt pair production, 


2.1. Observables and form factors 

In case of polarised beams Ref. [10] suggests to express the form factors introduced in Sec. 1 in 
terms of the helicity of the incoming electrons, 


= 


= 



( 4 ) 


with z = 1,2 and j = V,A and mz being the mass of the Z'^ boson. The tree level cross section for 
ti, quark pair production for an electron beam polarisation I = L, R reads 


O'/ 


2AN^P 


(1 + 0 . 57 -")(. FiV )" + (.F /;)2 + 
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where A = with the running electromagnetic coupling a(s) and is the number of quark 

colours. Furthermore 7 and /? are the Lorentz factor and the velocity of the t quark, respectively. The 
term describes the reduced sensitivity to axial vector couplings near the ti production 

threshold. The cross sections at the Born level of the signal process e+e“ —)■ ti and the main Standard 
Model background processes at a centre-of-mass energy of 500 GeV are summarised in Table 1. 


Channel 

^unpol. [ft>] 

0 --,+ [fb] 

(^+,- [fb] 

tt 

572 

1564 

724 


456 

969 

854 

u,d.s,c 

2208 

6032 

2793 

bb 

372 

1212 

276 


11185 

25500 

19126 

W+W- 

6603 

26000 

150 


422 

1106 

582 

Z^W+W- 

40 

151 

8.7 

^(T^O^D 

1.1 

3.2 

1.22 

Single t for e+e” —>■ e~vAb [11] 

3.1 

10.0 

1.7 


Table 1: Unpolarised cross-sections and cross-sections at tree level for 100% beam polarisation for 
signal and background processes. 

The forward-backward asymmetry can be expressed as 

{A%p)r = T AN,/3 ■ (g) 

O'/ 
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The sign applies in case of an initial left-handed polarised electron beam, i.e. I = L, and the ’+’ 
applies correspondingly in case of an initial right-handed polarised electron beam, i.e. I = R. In the 
Standard Model the forward-backward asymmetry takes the values {Apg)L = 0.37 and {App)ji = 0.45 
at tree level. 

Neglecting CP violating form factors, the fraction of right-handed t quarks is given by the following 
expression: 


{Fr)i — OAt :^{App)i (7) 

The sign applies in case of an initial left-handed polarised electron beam, i.e. I = L, and the ’+’ 
applies correspondingly in case of an initial right-handed polarised electron beam, i.e. I = R. The 
values expected in the Standard Model at tree level are {Fr)l = 0.25 and = 0.80. 

With the introduced observables the six CP conserving form factors defined for the Z'^ and the 
photon can in principle be extracted simultaneously. However, close to the tt threshold the observables 
depend always on the sum Fiy + F2V- Therefore, a full disentangling of the form factors will be 
imprecise for energies below about 1 TeV. Hence, in the present study either the precision on the 
Form Factors or equivalently on the Couplings gp p, are determined simultaneously, while the 

two F 2 V are kept at their Standard Model values or vice versa. Due to these considerations the study 
will only make use of the cross section and A^pp since these are either the most precise observable 
in case of the cross section or the one that is most sensitive to axial couplings in case of A^pp. It 
is however reminded that in [7] the fraction of right-handed t quarks is determined to a precision of 
about 2%. 

2.2. Theory uncertainties 

The extraction of form factors requires precise predictions of the inclusive top quark pair produc¬ 
tion rate and of several differential distributions. In this section the state-of-the-art calculations and 
estimate theoretical uncertainties are briefly reviewed. 

As discussed at the beginning of this section, the optimal experimental strategy should consider 
e“'"e“ —> W'^hW~b inclusively, without attempting to distinguish single top and top quark pair 
production. However, today, sufficiently precise calculations are not available for the full process 
e+e” —)■ W^hW~h. Therefore the discussion in this paper is based on the current state of the art 
calculations for e+e“ —)■ tt, assuming that in the next decade theorists will rise to the challenge of 
extending the calculations to e“''e“ —>■ W'^bW~b. 

The QCD corrections to e+e” — tt production are known up to N^LO for the inclusive cross 
section [12], and to NNLO for the forward backward asymmetry App [13]. The perturbative series 
shows good convergence. 

In the kinematic region at around ^/s = 500 GeV as relevant for this study the N^LO correction 
to the total cross-section is below 1 %. An estimate of the size of the next order - obtained from the 
conventional variation of the renormalisation scale by a factor two and one half - yields 0.3 %. It can 
therefore be concluded that the uncertainty of today’s state-of-the-art calculations is at the per mil 
level. 

In a similar manner the QCD corrections to the prediction of differential distributions and quantities 
such as the forward-backward asymmetry can be estimated. The size of the N^LO correction to A^pp 
is estimated using the scale variation to be smaller than 1%, see also the discussion in e.g. [13] . 

Electro-weak (EW) corrections to the same process have also been calculated. A full one-loop 
calculation is presented in [14] . The correction to the total cross section is found to be approximately 
5%. The electro-weak correction to the forward-backward asymmetry is large, approximately 10% [14, 
15] . Recent studies [16] show that the corrections are notably different for different beam polarisations. 
They change for example the shape of the angular distribution in case of 'Pf,-,'P^+ = —1,4-1 beam 
polarisation while they only influence the normalisation in case of Pf,-,Pf,+ = -1-1, —1. 

The above discussion refers to corrections to the process e“''e“ —)■ tt. Further corrections of order 
Tt/rrit ~ 1% are expected to appear if the decay of the top quarks is included in the calculation. 
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It can be concluded that the state-of-the-art calculations of QCD corrections offer the precision 
required for this study. Uncertainties are under relatively good control, with uncertainties to the cross 
section of the order of a few per mil and order 1% on the forward-backward asymmetry. Electro- 
weak (one-loop) corrections are large. Further work is needed to estimate the size of the two-loop 
correction and, ultimately, to calculate this contribution. Currently these aspects are discussed with 
theory groups. 


3. Analysis of simulated events 

The study has been carried out on a fully polarised sample albeit realistic values for the ILC are 
Ue-,Ue+ = ±0.8, =f0.3. The cross section and therefore its uncertainty scales with the polarisation in 
a well defined way according to [17] 

CT-P^_.-P^+ = ^ [(1 - ^e-7^e+)(o'-.+ ±cr+ _) ± (Ue- -T’e+)(cr+,- - cr_, + )] . (8) 

The observable varies only very mildly with the beam polarisation. The realistic beam polar¬ 
isation will be correctly taken into account in the uncertainty of the results. 

Signal and background events corresponding to a luminosity of 250 fb“^ at ^/s = 500 GeV for each of 
the two polarisation configurations are generated with version 1.95 of the WHIZARD event generator [18, 
19[ that provides lowest order calculations of the 2 —)■ 6 fermions subprocess and simulates multiple 
photon radiation from the initial state electron and positron in leading-logarithmic approximation. 
WHIZARD is interfaced to the PYTHIA Monte Carlo programme [20[ for QCD and QED showering. 
The generated events were subject to a full simulation of the ILD Detector [4[ and subsequent event 
reconstruction using the version ILD_ol_v05 of the ILD software. 

The decay of the t quarks proceeds predominantly through t —>■ W^h. The subsequent decays of 
the bosons to a charged lepton and a neutrino or a quark-anti-quark pair lead to a six-fermion 
final state. The study presented in this article focuses on the semi-leptonic final state l^vbhq'q. Several 
other Standard Model processes give rise to the same final state. The most important source is single 
t quark production. Another relevant source is the Z^W'^W~ production. Experimentally, Z^W^W~ 
production can be distinguished rather efficiently from t quark pair production. The separation between 
single t quark production and tt pair production is much more involved. Note however, that according 
to Table 1 single t quark production is strongly suppressed in case of 'Pf.-,'Pe+ = ±1,-1. 

The entire selection procedure including lepton and b jet identification, top quark reconstruction 
and suppression multi-peripheral 77 -^hadrons background is explained in detail in [7, 8[ and [21[. 
The total selection efficiency of about 56% for semi-leptonic tt events includes events with a r lepton 
in the final state. Background processes can be very efficiently removed down to a negligible level. 
A powerful tool is the b likeness or 6 -tag value that suppresses about 97% of the dominant W~^W~ 
background. Additional selection criteria comprise cuts on the t quark and boson masses and of 
the invariant mass of the total hadronic final state. 

With the determined efficiencies a statistical uncertainty of the cross section e"'"e“ —?► tt of 0.47% 
in case Ve-,'Pe+ = —0.8, ±0.3 and 0.63% in case Ve--,'Pe+ = ±0.8,—0.3 can be derived. 

3.1. Forward-Backward Asymmetry A*pg 

The forward-backward asymmetry A^pp has the well known definition 

,t _ NjcosOtop > 0 ) - NjcosOtop < 0 ) 

“ N{cos0top > 0) ± N{cos0top < 0) ’ ^ ^ 

where N is the number of events in the two detector hemispheres w.r.t. the polar angle 0top of the 
t quark calculated from the decay products in the hadronic decay branch. The direction measurement 
depends on the correct association of the 6 quarks to the jets of the hadronic W boson decays. The 
analysis is carried out separately for a left-handed polarised electron beam and for a right-handed 
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Figure 2: Left: Reconstructed forward-backward asymmetry compared with the prediction by the event 
generator WHIZARD [18] for two configurations of the beam polarisations. Right: The same but after 
the application of a on < 15 for the beam polarisations Ve-, 'Pe+ = ~lj +1 explained in the text. 
Note, that in both figures no correction is applied for the beam polarisations V ^-, ’Pe+ = +1) ~1- Tie 
figure on the right hand side shows also the residual Standard Model background. 


polarised beam. In case of a right-handed electron beam the direction of the t quark can be precisely 
reconstructed. In case of a left-handed electron beam the final state features two hard jets from the 
b quarks and soft jets from the hadronically decaying W boson. This configuration leads to migrations 
in the polar angle distribution of the t quark as visible in the left part of Fig 2. This implication 
motivates to restrict the determination of A*pg in case of 'Pe-,'Pe+ = ~1) +1 to cleanly reconstructed 
events. For this a test variable is defined that compares the measured values of the Lorentz factor 
7 of the top, the momentum of the b quark in the rest frame of the top and the angle cosObw between 
the b quark and the W boson. The reconstructed polar angle distribution of the t quark is compared 
with the generated one for different cuts on x^- For a value of x^ < 15 an excellent agreement between 
the generated and reconstructed polar angle distributions is obtained, see the right part of Fig. 2. The 
tight selection however reduces the efficiency in case of left-handed initial electron beams from 55% to 
28%. With this the forward backward asymmetry can be determined to a statistical precision of better 
than 2%. The precise results corrected to the beam polarisations Ve-,'Pe+ = ±0.8, =f 0.3 are given 
in Table 2 together with those for the cross section, see previous section. A more straightforward, 
albeit experimentally more challenging, way to control the migrations is to measure the charge of the 
b quarks that are issue of the t quark decay. References [22] and [23] describe the determination of the 
b quark charge using secondary tracks. The same value of Ap^ is obtained at a comparable selection 
efficiency [8[. This means that App can be determined with two independent methods. 


+ 

1 

{Sa/a)stat. [%[ 

{SApp/App)stat. [%[ 

-0.8,-P0.3 

0.47 

1.8 

-PO.8,-0.3 

0.63 

1.3 


Table 2: Statistical precisions expected for the cross sections and App for different beam polarisations. 

Hard gluon radiation may alter the polar angle distribution of the final state t quarks. The WHIZARD 
version 1.95 used for the study generates hard gluons only via the interface to PYTHIA that generates 
the parton shower. Therefore results presented before have been checked with a study on parton level 
using the most recent version 2.2.2 of WHIZARD that correctly accounts for hard gluon radiation. No 
significant difference has been observed. 
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4. Discussion of systematic uncertainties 

In the previous sections measurements of either cross sections or asymmetries have been presented. 
This section makes an attempt to identify and quantify systematic uncertainties, which may influence 
the precision measurements. 

• Luminosity: The luminosity is a critical parameter for cross section measurements only. The 
luminosity can be controlled to 0.1% [24]. 

• Polarisation: The polarisation is a critical parameter for all analyses. It enters directly the 
cross section measurements. The studies presented in [25] using W pair production lead to an 
uncertainty of 0.1% for the polarisation of the electron beam and to an uncertainty of 0.35% for 
the polarisation of the positron beam. This translates into an uncertainty of 0.25% on the cross 
section for Ve- , Ve+ = —0.8, +0.3 and 0.18% on the cross section for Ve- , 'Pe+ = +0.8, —0.3. The 
uncertainty on the polarisation can be neglected with respect to the statistical uncertainty for 

^FB- 

• Beamstrahlung and beam energy spread: The mutual influence of the electromagnetic fields of 
the colliding bunches provokes radiation of photons known as Beamstrahlung. This Beam¬ 
strahlung modulates the luminosity spectrum, i.e. moves particles from the nominal energy to 
smaller energies. At the ILC for a centre-of-mass energy of 500 GeV about 60% of the particles 
are expected to have 99% or more of the nominal energy [4[. The beam energy spread, i.e. the 
RMS of this main luminosity peak is 124 MeV for the electron beam and 70 MeV for the positron 
beam [4[. Both effects play a role at the tt threshold [26] and can be neglected at energies well 
above this threshold. 

• Experimental uncertainties in top quark reconstruction: As discussed in Sec. 3.1 migrations have 

to be taken into account for the measurement of in particular for the polarisations Vg -, Ve+ = 
—0.8, +0.3. These migrations are reduced by stringent requirements on the event selection using 
a analysis. This in turn leads to a penalty in the efficiency. The success of the method depends 
in addition on a very sharp measurement of the variables used for the analysis. It is expected 
that these ambiguities can be (partially) eliminated by an event-by-event determination of the 
charge of the b quark from the t decay. As has been shown in Sec. 3, the effect will be very much 
suppressed in case of = +0.8, —0.3 beam polarisation. 

• Other experimental effects: There is a number of other experimental effects such as acceptance, 
uncertainties of the b tagging or the influence of passive detector material. The LEPl experiments 
quote a systematic uncertainty on Ri, of 0.2% a value which may serve as a guide line for values to 
be expected at the ILC, which on the other hand will benefit from far superior detector resolution 
and b tagging capabilities. 

• Theory: The uncertainties of today’s state-of-the-art calculations are discussed in Section 2.2. 
The uncertainties in the QCD corrections to the total cross section and Afb are of the same order 
as the experimental uncertainties. Two-loop electro-weak calculations are required for a reliable 
estimation of the uncertainties due to electro-weak corrections. It is however intuitively clear that 
the latter will benefit from the insight of the different impact for different beam polarisations, 
see Sec. 2.2 and [16[. 

• Single-top production: Single top production at the LC in association with a W boson and bottom 
quark (through WW* production) leads to the same final state as t quark pair production. 
Being largely suppressed in case of 'Pe-,'Pe+ = +0.8,—0.3 beam polarisation, it forms a sizeable 
contribution to the six-fermion final state in case of Ve-, Ve+ = —0.8, +0.3 beam polarisation. It 
must therefore be taken into account in a realistic experimental strategy. This is left for a future 
study. 
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• Beyond Standard Model Physics: Possible BSM effects may affect the various components of 
the background, in particular the tt induced background. This will therefore require a careful 
iterative procedure with tuning of event generators. This procedure seems feasible without a 
significant loss of accuracy. 

As a summary it can be concluded that the total systematic uncertainties will not exceed the 
statistical uncertainties. This, however, requires an excellent control of a number of experimental and 
theoretical quantities. 


5. Precision of form factors and electro-weak couplings 

The measured cross sections and lead for two polarisation configurations to a set of four 

observables. By means of Eqs. 5 and 6 the uncertainties on these observables are used to build up 
a system of linear equations to determine the variances of up to four variables^, The variances are 
equivalent to the square of the standard deviations of the variables under study. These variables can 
be the form factors or alternatively directly the couplings. More explicitly, in this paper the following 
quantities will be determined separately: 

1. The standard deviations of the Form Factors F^y, F^y, F^^ assuming no variation of the Form 
Factors F^- 

2. The standard deviations of the Form Factors F^y, F^y assuming no variation of the Form Factors 

px . 

3. The standard deviations of the Couplings 

Note, that the Form Factor Ffy^ is fixed to be 0 in order to respect QFD gauge invariance. On the 
other hand all four Couplings gf a,ie allowed to vary freely. The resulting standard deviations are listed 
in Table 3. 


Quantity 


F'^ 

F^ 

^2V 

F^ 

9l. 

9r 

9f, 

9r 

SM Value at tree level 

2/3 

0.230 

-0.595 

0 

0 

2/3 

2/3 

0.824 

-0.364 

Standard deviation 

0.002 

0.003 

0.007 

0.001 

0.002 

0.005 

0.005 

0.008 

0.009 

Relative precision [%\ 

0.3 

0.9 

1.2 

- 

- 

0.8 

0.8 

1.0 

2.5 


Table 3: Standard deviations and resulting relative precisions of form factors and couplings derived 
from the statistical precisions on the observables cross section and Ap^ as listed in Table 2. 

The complete covariance matrices are given in Appendix A. From there it can be told that e.g. gf and 
gp are strongly correlated while Ffy and are nearly uncorrelated. The expected high precision 
at a linear e“''e“ collider allows for a profound discussion of effects by new physics. The findings can 
be confronted with predictions in the framework of Randall-Sundrum models and/or compositeness 
models such as [2, 34, 36, 32, 31, 28, 29, 33] or Little Higgs models as e.g. [30]. All these models 
entail deviations from the Standard Model values of the t quark couplings to the boson that will 
be measurable at the ILC as illustrated in Fig. 3. Therefore, the couplings of the t quark to the 
boson will discussed in a broader context in the following. 

5.1. Discussion of couplings to the Z^ Boson - Comparison with perspectives for LHC and Flavour 
Physics 

Flectro-weak couplings can be measured at the LHC in associated tty and ttZ production. A 
comprehensive compilation on the statistical precisions on the form factors that can be expected at 
the end of the HL-LHC is given in [37] and [38] for an update on ttZ^ form factors. The values 



















20 % 


10 % 


RS with Z-Z’ Mixing [34] 


Light top partners Alternative 2 [29] 



SM 


-330% 


- 20 % 


- 10 % 


Light top partners 
Alternative 1 [29] • 

5D Emergent [32] 


- 10 % 


- 20 % 


+ 


+ 




10 % 20 % 

• 4D Composite Higgs Models [33] 

Light top partners [28] 

Little Higgs [30] 

RS with Custodial SU(2) [31] 

Composite Top [2] 


Figure 3: Predictions of several models that incorporate Randall-Sundrum (RS) models and/or com¬ 
positeness or Little Higgs models on the deviations of the left- and right-handed couplings of the 
t quark to the boson. The ellipse in the frame in the upper right corner indicates the precision 
that can be expected for the ILC running at a centre-of-mass energy of ^/s = 500 GeV after having 
accumulated L = 500 of integrated luminosity shared equally between the beam polarisations 
Vg -, Ve+ = ±0.8, =f0.3. The original version of this figure can be found in [35]. 
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Figure 4: Graphical comparison of statistical precisions on CP conserving form factors expected at 
the LHC, taken from [37] and [38], and at the ILC. The LHC results assume an integrated luminosity 
of C = 3000 fb~^ at y/s = 14 TeV. The results for the ILC assume an integrated luminosity of 
C = 500 fb~^ at y/s = 500 GeV and a beam polarisation Ve- = ±0.8,7^6+ = tO.3. 


published there are compared with the results in the present study in Fig. 4. All but one form factor 
will be measured at about a factor 10 better at the ILC for the scenario discussed in this paper than 
it will be possible at the LHC. This exception is where [38] quotes a possible statistical precision 
of « 0.031. It should however be pointed out that the considerable precision expected for 6F^j^ 
benefits strongly from LEP/SLC bounds on the oblique parameters that e.g. render it unlikely that 
flips sign due to New Physics. The study presented by [38] is an analysis at leading order QCD. 
The analysis carried out in [39] suggests that higher-order effects in the theory may allow for an 
improvement of the LHC precision by up to 40%. Note at this point that the interference between the 
7 and the in case of e+e“ —?► ti will allow for measuring flips of the signs of the form factors that 
will be unnoticed in associated ttZ^ at the LHC. 

While the prospects for the LHC discussed so far are based on analyses differential in given jet 
observables of the final state, LHC experiments observe the process pp itZ^ [40, 41, 42, 43] . The 
interpretation of the results is however still limited by the small statistics available for the analyses. 

At the LHC electro-weak couplings are measured also in single t quark production. In the effective 
field theory approach, assuming SU{2)i^ x U{1)y gauge symmetry for the operators, the relation 


„tbW 

SL 


0.35 




( 10 ) 


can be established. Here is the charged current coupling of the decay t —>■ Wb. The CMS 
Collaboration [44] reports a precision for the t-b transition probability Vtb of about 4%. In the Standard 
Model Vtb is identical to g*^^. Hence, by means of Eq. 10 the precision of the coupling of left-handed 
t quarks to the Z boson can be derived to be of the order of 11%. Noting that <j{pp —>■ itZ) ^ 


Wor the Linear Algebra the software package Eigen [27] version 3.2.2 has been used. 
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+ (Sr)^ allows in principle also for deriving ( 5 ^)^, albeit with a poor precision given that 
^ Loop corrections in heavy flavour physics as e.g. in the processes b —>■ sy, B 

01 K ^ fX^, respectively, may also lead to competitive determinations of Sgf^ [45]. However, again 

can only be constrained rather poorly. 

It follows that the ILC will allow for measurements superior to those that can be expected from 
existing experiments. This is particularly true for the determination of 6g^. 

5.2. Coupling measurements and form faetors at different centre-of-mass energies 

Given the fact that at the ILC in its current layout centre-of-mass energies of up to 1 TeV can 
be reached and that the alternative project for a linear collider, CLIC [46], may even reach higher 
energies, it is instructive to discuss the results presented in this paper with this possibility in mind. 
The selection and reconstruction of the decay topology of boosted t quarks is very different from that 
of t quarks with moderate velocity. Therefore, the study must be extrapolated to high centre-of-mass 
energy with some care. Still, the following observations can be made: 

• Neglecting varying detector systematics and theory uncertainties with varying centre-of-mass en¬ 
ergy, and assuming the linear collider luminosity vs. centre-of-mass energy curve, the sensitivity 
to the form factors considered in this paper is greatest at approximately ff~s = 400 — 700 GeV. 
At lower centre-of-mass energy, as e.g. studied in [47], the small velocity of the t quarks reduces 
the potential of the measurement compromising thus the measurement of the axial vector 
couplings to the boson and by virtue of Eq. 3 the disentangling of left- and right-handed 
couplings. On the other hand running at centre-of-mass energies close to the tt threshold offers 
sensitivity to virtual Higgs exchange [48, 49, 50] . In case the Higgs has a CP odd component this 
may give rise to recognisable CP violating effects in the threshold region [51]. However, in the 
transition region between the ti threshold and the continuum region starting at around 380 GeV 
the current QGD uncertainties are at least 10%. This is due to uncertainties on higher QGD 
order corrections and on the correct matching procedure between the non-relativistic calculations 
at the ti threshold and the relativistic continuum calculations [52] . 

• If an effect is seen at ^/s = 500 GeV it will be crucial to know how it evolves with energy with a 
decent lever arm. If, for instance, the effect is due to mixing of the boson with a new Z' boson 
it will remain unchanged. If, however, a Z' boson leads to a propagator term, the corresponding 
effect will grow like s/M"^,. In the case of Randall-Sundrum Models both effects are present 
and therefore measurements at two energies are needed to extract Mz', see e.g. [35] for a deeper 
discussion. 

• The impact of high-scale new physics on the observables can increase strongly with centre- 
of-mass energy. Operators corresponding to the top quark dipole moments and four-fermion 
contact interactions induce larger anomalous form factors at higher energy. For other anomalous 
couplings, however, the impact is nearly independent of the centre-of-mass energy as is the case 
for F^y and 

A full simulation study at different centre-of-mass energies is left for a future publication. 

6. Summary and outlook 

This article presents a comprehensive analysis fully simulated events of ti quark production at the 
International Linear Gollider using the semi-leptonic decay channel. Results are given for a centre-of- 
mass energy of ffs = 500 GeV and an integrated luminosity of £ = 500 fb“^ shared equally between 
the beam polarisations Ve-, Pe+ = ±0.8, =f0.3. 

Semi-leptonic events, including those with t leptons in the final state can be selected with an 
efficiency of about 55%. The cross section of the semi-leptonic channel of ti quark production can 


11 


therefore be measured to a statistical precision of about 0.5%. The second observable is the forward- 
backward asymmetry . It was shown that in particular for predominantly left-handed polarisation 
of the initial electron beam the V — A structure leads to migrations, which distort the theoretical 
expected A^pg. These migrations can be remedied by tightening the selection criteria of the events or 
alternatively by measuring the charge of the b quark produced in the decay of the t quark. Taking into 
account this correction the forward-backward asymmetry can be determined to a statistical precision 
of better than 2% for both beam polarisations. 

The observables together with the unique feature of the ILC to provide polarised beams allow for 
a largely unbiased disentangling of the individual couplings of the t quark to the boson and the 
photon. These couplings can be measured with high precision at the ILC and, when referring to the 
results in [37, 38], considerably better than it will be possible at the LHC even with an integrated 
luminosity of £ = 3000 fb“^. The improving analyses of the LHC experiments will however be observed 
with great interest. 

Beam polarisation is a critical asset for the high precision measurements of the electroweak t 
quark couplings. Experimental and theoretical effects manifest themselves differently for different 
beam polarisations. It seems to be that the configuration Ve-, 'Pe+ = -1-0.8,—0.3 is more benign in 
both, experimental aspects due to the suppression of migrations in the polar angle spectrum of the final 
state t quark and theoretical aspects due to the somewhat simpler structure of higher order electroweak 
corrections. It is intuitively clear that the described facts would greatly support the discovery of effects 
due to new physics. 

The precision as obtained in the present study for the ILC would allow for the verification of 
a great number of models for physics beyond the Standard Model. Examples for these models are 
extra dimensions and compositeness. The results obtained here constitute therefore a perfect basis for 
discussions with theoretical groups. Note at this point that running scenarios for the ILC have been 
proposed that would yield between 8 and 10 times more integrated luminosity [53] than it is assumed 
for the present study. Moreover it can be expected that the event reconstruction will be improved 
by e.g. the measurement of the b quark charge. It is therefore not statistics that will limit the final 
accuracy but most likely theory and experimental systematics. 

Hence, the study of systematic errors, only partially addressed in this study, will become very 
important. Already from the achieved precision it is mandatory that systematics are controlled to the 
1% level or better in particular for the measurement of the cross section. This issue is addressed in 
ongoing studies. 

The study presented in [16] based on generated events suggests that by exploiting the polarisation 
of the final state t quarks a simultaneous extraction of all ten form factors, see Eq. 1, to a precision 
below the percent level is feasible. A detailed comparison between the advantages and drawbacks of 
the method applied there and the method presented in this paper is left for a future study. 
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A. Covariance matrices 

For completeness the underlying covariance matrices of the results presented in Sec. 5 are given in 
this appendix. 
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• The covariance matrix resulting from the system of linear equations built for the Form Factors 
PiA reads: 


var{F2y) cov{F^y,Ffy) cov{F^y, Ffj^) 


■ 0.260 

- 0.043 

0.506 ■ 

var{F^y) cov{F^y,Ff^) 

= 


0.791 

0.118 

var{Ff^) 




5.460 


• The covariance matrix resulting from the system of linear equations built for the Form Factors 
F^y reads: 


var{F]y ) cov{F]y, F^y ) 
var{F§y) 


0.160 - 0.046 
0.473 


X 10"^ 


( 12 ) 


• The covariance matrix resulting from the system of linear equations built for the Couplings 
91.9r,91.9r reads: 


var{gl) 


cov{gl,gJf) cov{gl,gl) 
var{gl) cov{g},,gl) 
varlgl) 


2.801 - 2.357 - 0.339 
2.949 - 0.780 

5.988 


cov{gl,g^) 

cov{g]^,g^) 

covigl,g^) 

var{g^) 


- 1.234 

2.247 

- 4.691 


X 10"^ 


7.306 


( 13 ) 
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